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Disclaimer 



/ 

Commercial materials are identified in this report in order to 
adequately specify the experimental procedure. In no case does such 
identification imply recommendations or endorsement by the National 
Bureau of Standards, nor does it imply that the equipment identilied 

is necessarily the best available for the purpose. j 

j 
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Abstract 


The development of an experimental small-scale modeling technique 
for conducting simulations of air distribution patterns in full-sized 
rooms is described. The technique was applied to simulate the motion 
of air in a single-room survival shelter. The parameters that describe 
and must satisfy the principles of dynamic similarity for a comparison 
of convection between a model and a full-scale room were determined. 

The main parameters were the Grashof and Prandtl numbers of the fluid 
and the aspect ratio of the room. High temperature water was used in 
place of air as the fluid in the model. The model was 6" x 6" (width 
oy height) and 18" deep. A tracer particle was injected into the 
water in the model. Time exposure photographs were made of the tracer 
to reveal fluid motion. A description of the experimental apparatus 
is given including a sample test procedure. Experimental results ob¬ 
tained are presented and discussed. 
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1. Introduction 


In an emergency, survival shelters may be occupied for several dav 
or perhaps for several weeks. The actual occupancy period must depone 
on conditions prevailing outside the shelter. The comfort, indeed the 
survival, of an occupant in such a shelter depends on whether 01 not 
the shelter receives an adequate amount of ventilation and if the 
ventilation air is properly distributed within the shelter. Since 
electric power cannot be depended upon completely for mechanical venti¬ 
lation during the time when the shelters are likely to be occupied, 
the forces that produce natural ventilation should be used to a maximum 
Some research has been reported on natural ventilation of both 
above and below ground shelters [1, 2], This work concentrated mainly 
on the quantity of air moving into and out of a shelter with little 
emphasis on the distribution of air within the shelter. Studying the 
air distribution in a full scale shelter would be both time consuming 
and costly; particularly when a variety of conditions are considered. 
For this reason, the major emphasis of this program was given to the 
development of a small scale modeling technique which would be applica¬ 
ble for naturally ventilated survival shelters. 





The modeling of any convection system requires close agreement be¬ 
tween the similarity parameters of the model and tlu- actual system being 
simulated. For a natural convection system these parameters are usually 
limited to Grashof and Prandtl numbers ami a geometrc dimension. The 
premise upon which this study is based is that because of the agreement 
between their Grashof and Prandtl numbers and aspect ratio high tempera¬ 
ture, high pressure water in a small "table top" mod.'I adequately simulates 
convection patLcrns of air in a full size! room. 

The objecL vc of t: s study was ultimately lo obt> overall con¬ 
vection flow descriptions for single room shelters for both vented and 
non-vented conditions, with and without occupants. However, since the 
experimental method offered a number of advantages for future work as 
well, the development of the method was an equally important objective. 

This report is presented in two distinct parts. The main body considers 
the theoretical justification of the model and lie actual experimental 
results. The Appendix is a translation of these results to a simulated 
case of air in a full sized room. Although this translation has several 
limitations it Is felt that the results yielded arc of significance to 
the designer of survival shelters. 


2 




2. Theoretical Considerations for Scale Modeling 

The differential equations governing continuity and the transfer 
of momentum and energy are well known for the case of natural convection 
in a space enclosed by a rectangular parallelepiped. A study of these 
equations was made in order to determine modeling criteria. 

The general differential equations and boundary conditions governing 
the behavior of a fluid within an enclosure may be oppressed as folLows 
[5] : 

a. Momentum equation 

DW 

p — = F - grad p + p ; r W 


b. Energy equation 


pC 


91 = 

Dt 


div W . 
P “J 


- div <; 


c. Continuity equation 


m <e w > * 0 

I 

d. Boundary conditions along the inside surfaces of the 
enclosure may be described as 

W = 0 

T; a given function of position 
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In the above equations the following nomenclature is employed: 

W = velocity vector 
T = temperature 
P = pressure 

F = external force (such as gravitational force) vector 

.. = dynamic viscosity 

o = density 

C = specific heat 

* = thermal condo. *ivity 

J = mechanical equivalent of heat 

*' = Rayleigh energy dissipation function 

q = radiation heat flux vector 
r 

t = time 

In the energy equation the term involving ,> div W is the rev tsihU 
compression work done by the fluid, and the term involving a © i: the 
irreversible work done on the fluid element. These terms a: usually 
very small for the natural convection process. 

Assuming that the convection system can be characterized by a 
reference length d (such as the width of the enclosure) and a reference 
temperature difference ' (the difference between the maximum surface 
temperature, T, , and minimum surface temperature, T ), the following 
dimensionless variables may be introduced: 
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dimensionless velocity, U = — 

where kinematic viscosity •> = i /,• 


dimensionless temperature, H 


T - T. 


dimensionless pressure, P = - 


dimensionless time, t - Jr 


dimensionless external force, f - F/(p - r ) g 


dimensionless length, X - x/d 

Y = y/d 


The differential equation can now be expressed with these dimen¬ 
sionless parameters as follows: 


^ = Gr • 0 * f - grad P + " U 


D0 _ 1 , P _ 1_ 

Dt Pr pC M r 


where Gr and Pr are Grashof and Prandtl numbers respectively and defined 


as follows 


Gr = 


and 0 = thermal expansion coefficient of the convecting fluid. 
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A 


The dimensionless forms of the differential equations suggest that 
a solution to this set of differential equations can he applicable to a 
variety of natural convection systems of different thermophysical proper 
ties and characteristics. In order for this to be so, each system must 
have the same values of Gr and Pr, have the same radiation vei l.wr 
characteristics, and have similar geometrical shape and temperature 
profile along the surface boundaries. 

For example, consider two convection systems with characteiist. ic 
lengths, d-j and d , reference temperature di 11 eroncos - and 6 , and 
the reference temperature levels 1’ » and T , . t ho velocity and ' •■mpera- 

ture patterns of both systems are correlated is follow;: 


II = 



l-’ d. 



1 


for time t. and t , such that 

v. t. . t r 

dw d. 

The same relationships may be used to translate experimentally observed 
velocity and terperature patterns of system 1 into those for system 2. 
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V 

4 The Grashof number for natural convection in a room 8' wide by 8' 
high and with a 10 ®F temperature difference between opposite walls is 
on the order of 10^ (based on the width, d, as the critical dimension). 
This is within the region of turbulent convection since the regior of 

9 

transition from laminar to turbulent fiow is on the order of (10 ). The 
Prandtl number of air at 80 °F is about 0.72 |6], and the aspect ratio 
of the above mentioned room is 1. 

In order to model the natural convection of air in a room, it is 
desirable to have a fluid which will produce the same Grashof number and 
Prandtl number - in a scale model as air will have in an average-size 
room. A number of fluids were found which would give a high Grashof 
number in a small (6" high x 6" vide) enclosure, but their Prandtl number 
was much higher than that for air. Of the many fluids considered, only 
high temperature water gives a high Grashof number and a Prandtl number 
close to that of air (see figure 1). Water at 500 °F has a Prandtl 
number of 0.87 [6] which approaches the value of .72 for 80 °F air. 

Figures 2, 3, and 4 show the thermophysical properties of air and 
water that comprise the Prandtl number. The Prandtl number of water 
decreases rapidly from a value of 6.5 at room temperature to slightly 
less than unity as it is hefted to 500 °F. This reduction is primarily 
due to a decrease in viscosity. Examination of these figures provide 
the following important information: 

a. High temperature water, at 500 °F, used inside a 
small-scale model can provide a convection fluid 
with the Prandtl number close to that for convec¬ 
tion conditions found in room air. 
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b. Thertnophysical properties of hot water are relatively 
constant with respect to temperature change, thus 
inaccuracies due to property value variation across 
the boundary layer along the heat transfer surfaces 
are minimized. In contrast the viscosity of low tem¬ 
perature water changes considerably with the tem¬ 
perature changes across a boundary layer, severely 
affecting the convection flow pattern. 

Although not shown in these figures, the pressure level for the 
high temperature water tank increases very rapidly with temperature. 

Thus high pressure and high temperature must both be considered when 
designing an apparatus in which high temperature water is used as a 
modeling fluid. 

It was concluded that natural convection in a room can be simulated 
by a small-scale model using hot water in place of air. This conclusion 
is valid if the effect of thermal radiation can be ignored in the original 
differential equations, or if the radiation effect is equal for hot water 
and air. A detailed analysis of this subject is involved and a compre¬ 
hensive treatment can be found in reference [5]. The assumption that 
convecting fluids are gray (radiant heat transfer independent of wave¬ 
length) and non-scattering will be discussed here briefly. 
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In the absence of carbon dioxide and water vapor the air in a room 
may be considered transparent to thermal radiation. But carbon dioxide 
(up to 37.) and water vapor (up to 1.57. at 70 °I ) are pre. ent in the at¬ 
mosphere and they do absorb some of the incident radiant energy. Also 
airborne dust particles scatter some of this incident radiation. For ex¬ 
ample, the absorbing and scattering effects of the air surrounding the 
earth reduce by 207. the amount of normal incident solar radiation at the 
earth's surface. This holds true even during cloudless days (the ap¬ 
proximate absorption coefficient for air is 0.374 ft ^ at standard at¬ 
mospheric pressure). 

However, since the optical thickness involved ip a room size en¬ 
closure is very much smaller than the atmospheric layer, it has been 
customary to ignore the radiation term fdiv q^) in .he energy equation. 
The rationale for this is the fact that the experimentally observed air 
velocity and temperature profiles in the laminar boundary layer flowing 
along a vertical plate show good agreement > ith the calculated values 
of velocity and temperature based upon the energy equation without the 
radiation term [7]. 

On the other hand, liquid water is practically opaque to thermal 
radiation (its absorption coefficient is approximately 15,000 ft ^). 
According to Sparrow [5], this constitutes a case for the optically thick 
limit, which in effect permits the div q^ term to be omitted. 
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Neglecting the radiation term in the energy equation, however, is 
not meant to imply that radiant heat transfer is unimportant in room air 
convection. On the contrary, room air motion is profoundly affected by 
the radiant heat exchange among room surfaces by healers and by incident 
solar and sky radiation entering through windows. However, the influence 
of this radiant heat transfer first appears in the form of a surface tem¬ 
perature pattern in the room rather than as a direct influence or some 
fluid element within the room. It is this surface Lemperaturc pattern 
then which actually influences convective heating and the subsequent 
change of air motion. This means that the convcctiv ■ '>oat transfer from 
room surfaces to room air and the radiation heat exchange among various 
surfaces must be considered simultaneously. 

Such a combined analysis of heat transfer mode.. extremely com¬ 
plex and probably iterative hat radia it heat transfer among surfaces 
in a room can be calculated independently. The diilinulty caused by 
radiant heat exchange among room surfaces can, however, be circumvented 
experimentally by insuring that surface temperatures of a small-scale 
model be kept very close to those of a room being modeled. 

Although high temperature water is an excellent fluid for the scale- 
model study of room air convection, the high pressure associated with this 
fluid becomes a major consideration in the design of an experimental ap¬ 
paratus. A low Prandtl number of water f0.87) can be attained only by 
heating the water to 500 °F requiring a saturation pressure of 680 psia 
to prevent boiling. The tank must safely restrain high temperatures and 
pressures, provide an adequate size view window for flow observations, 
contain a system for both measuring and controlling temperature, and 
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accommodate a mechanism for injecting the trace particles. 

There is yet some question as to whether or not the flow and tem¬ 
perature patterns are truly dependent upon the Prandtl number when turbu¬ 
lent flow conditions actually occur in a room or shelter. Moreover, it 
is hypothesized that the heat transfer rate at the fluid-surface inter¬ 
face is less affected by the Prandtl number value for turbulent boundary 
layer flow than it is for laminar flow condition. 

To evaluate the magnitude of the effect of Pr upon convection 
characteristics, the energy equation (without the radiation term) was 
solved numerically with the use of a finite-difference technique [8]. 
Figures 5, 6, 7, 8, 9, 10, 11 and 12 illustrate some of the results ob¬ 
tained. As showi, Pr and Gr strongly influence both velocity and tern- 

3 6 

perature profiles for a span Gr from 10 to 10 . The convection flow 

at both these Grashof numbers are in the laminar region. Numerical 

£ 

calculations for Gr higher than 10 are extremely time- and memory¬ 
consuming for a computer, but they are not impossible. 

Rockett and Torrance [9, 10] reported a numerical and experimental 
study of convection of air in a room. In their work, however, air 
movement was caused by a localized source of heating, such as a waste¬ 
basket fire. They used a "finite-difference" technique to solve the 
equations governing natural convection. They also conducted an experi¬ 
mental study using air in a scale model and found close agreement between 
the numerical and experimental results in the laminar-flow range. This 
close agreement adds confidence to the theory that the Grashof and Prandtl 
numbers of the fluid and the aspect ratio of the enclosure are the factors 
which must be kept consistent between a scale model and a full-sc^le room 
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in order to show dynamic similarity of the flows in each. 

The Rockett and Torrance studies were carried out mainly in the 
laminar flow range. Only one of their numerical simulations was made in 
the turbulent range and it is difficult to judge how close an agreement 
was obtained between experiment and theory. Thus, the question remains 
as to the real necessity of having similar Pr and Gr numbers between 
small scale model and full sized room, in the turbulent range of con¬ 
vection. 

An approximate analysis for the turbulent region can be made by a 
von-KaVman-type solution as was done by Eckert and Jackson [11]. They 
analyzed the experimental data of Griffiths and Davis [12] and found 
that the velocity and temperature profile of the turbulent natural con¬ 
vection boundary layer along a heated vertical flat plate can be ex¬ 
pressed by the following equations: 

1/7 4 

" ’ ", <f> 

e « 9 (i - J) 

where 

U x = 1.185 £ (Gr) 1/2 <1 + O.Oi94 (Pr) 2/3 )" 1/2 

6 = 0.565 x (Gr)" 0,1 (Pr)" 8/l5 (1 + 0.0494 (Pr) 2/3 )" 1/2 

y is the coordinate perpendicular to the plate 
x is the coordinate parallel to the plate 
® is the temperature difference between the wall and the 
free stream 

6 is the boundary layer thickness at x. 
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In order to examine the combined effect of U T and 6» the net total flow 
may be calculated as follows: 


6 1/7 4 

P J {1 " dy 

o 

Making use of the dimensionless notations introduced previously a dimen¬ 
sionless total flow, q, can be defined as 


q = f UdY = C, (Gr) 


0,4 (Pr)~ 8/15 (1 + 0.0494 p r 2/3 )“ 0 * 4 


o 

where C, is a constant and its exact value is not important for this 
discussion. 

A similar relation for the laminar boundary layer flow or for the 
convection with the low Grashof number is 


q = J UdY = C a (1 + 1.05 Pr)“° -25 (Pr)'° -5 (Gr)°' 25 

A calculation was then made to check the influence of the Prandtl 
number on the flow in a rectangular cavity. When these net flow parameters 
were evaluated at the laminar Grashof number, the overall flow rate for 
Pr = 0.7 was found to be 4.44 times as much as that for Pr = 7.0. While 
at a turbulent flow condition the overall flow rate for Pr = 0.7 was 
3.55 times as much as for Pr = 7.0. 

This analysis, therefore, ind*cr.twS that the Prandtl number dependence 


on the convection flow pattern is still very strong in a condition of 
turbulent boundary layer flow and high Grashof number. 





3. Experimental Apparatus 


A steel pressure chamber was constructed, to contain high temperature 
water, as shown schematically in figure 13. The overall dimensions were 
20" long by 12" in diameter. The viewing window was of tempered plate 
glass 3/4" thick and 8" in diameter. The lighting port at top was 1" 
wide and 10" long. 

Figure 14 is a schematic of the modeling enclosure which was located 
inside the chamber. It was a rectangular parallelepiped 6" x 6" x 18" 
with the front end open for observing the flow. It was made of 1/8" 
thick copper faced with black anodized aluminum plates to provide good 
contrast in photographs. A slit was cut in the ceiling of the model, 
in line with the lighting port of the chamber, and fitted with a glass 
for passing light to the interior. On the floor of the enclosure, also 
in line with the lighting port in the outer chamber, was located a mirror 
to increase illumination. On the outside of the enclosure serpentine 
coils were soldered to the walls and floor for circulating water to 
allow individual control of surface temperature. 
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Figure 15 Schematic of Modeling Enclosure Inside the Pressure Chamber 











In figure 14 thermocouples are shown behind the light slot, that is, 


toward the rear of the model farthest from the viewing window. Eight 
thermocouples were held against the surfaces, two on each of the four 
surfaces. Sixteen of them were held away from the walls, positioned in 
a grid to treasure the fluid temperature, and two were used to measure 
the temperature of water in the two ventilation ducts. The temperatures 
were read by an indicating potentiometer. The thermocouple leads were 
40-gage copper-constantan wire jacketed in a stainless steel tube. They 
entered the chamber through a pressure fitting in the removable rear 
plate. 

Figure 15 shows the enclosure as it fit inside the pressure chamber. 
Figure 16 is a schematic showing the heat exchange system used to provide 
hot water to the serpentine coils and to the natural ventilation heat 
exchanger. The modeling enclosure was ventilated naturally by a supply- 
and-return duct. These ducts were connected to the natural ventilation 
heat exchanger and entered the top of the modeling enclosure through the 
light entry port. They were used to simulate the exchange of fluid with 
an external environment. The portion of both ducts directly above the 
light entry port was made of glass tubing to minimize shadows in the 
modeling area. 














A rust inhibitor, sodium chromate, was added to the water in the 
chamber. About 0.8 oz. was used normally. The chromate added a yellow¬ 
ish cast to the water, but this was not objectionable since it did not 
seem to impair the quality of the photographs until the temperature rose 
above 300 °F. Above 300 °F the water became a rather reddish brown color 
and was too cloudy for good photographs. (Sodium chromate should be 
handled with caution because it is quite toxic.) 

The tracer particles were small aluminum flakes 4 x 3 x 10 p) 
which had been sifted through a 325 mesh during manufacture. A coating 
of styric acid for lubrication was applied to the flakes. The tracer 
injector is shown in the schematic in figure 17. This device was needed 
because the tracer particles had to be injected at a pressure slightly 
above that of the chamber. The injector was designed so that the pres¬ 
surizing gas was not carried into the chamber with the tracer since 
this would cause bubbles which would impair the photographic results. 

The simulated occupant, which can be seen in the photographs toward 
the end of the Results Section, contained a heating element which was 
used to maintain its surface temperature above that of the surrounding 
water. 
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Schematic of Lighting 



The lighting equipment is also shown in the schematic on figure 17. 
The lights were high-intensity photo flood lights. Variable voltage 
transformers controlled the power to the lights. Ammeters were placed 
in the lighting circuits so that identical intensity levels could be 
repeated during photographic sessions. The collimating slits shown 
between the lights and the light entr/ port on the chamber were ad¬ 
justable in width and in vertical position. A 1/2" slit width usually 
provided adequate illumination in the model. 

The camera shown in tlgure i3 was mounted on a base having three 
degrees of freedom for proper position and focusing. The camera had an 
f2.8, 100 nm lens. The back was detachable so that a ground glass fo¬ 
cusing screen, a roll film back, or a Polaroid back could be used. The 
roll film used was 120 Trl-X-Pan film. 

4. Experimental Procedure 

As a first step the entire chamber was flushed with water. The 
modeling area, especially, was cleaned of aluminum particles used during 
the previous test. The flushing removed not only the trace particles 
but also the cloudiness due to rusting which was present after each 
test. Often three or four flushings were required to thoroughly clean 
the chamber The viewing window, the light entry slot at the top of 
the modeling area, and the light entry port at the top of the chamber 
were cleaned with an ordinary glass cleaning compound. 
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The apparatus was designed so that the majority of the necessary 
changes within the model could be made through the front viewing port 
which was easily removed. Most of these changes involved varying the 
positions of the ventilation ducts and installing simulated occupants. 

After the last flushing the chamber was sealed, and hot water added 
through the fill line. Hot water was used because less time was then 
needed to reach normal operating conditions, and the hot water usually 
contained less dissolved air. Sodium chromate was added to this water. 

Once the chamber was filled, the pump and the heaters were turned 
on to bring the water in the chamber to operating temperature. By the 
time the water in the model reacnea 205 “F, or slightly below boiling, 
the air originally dissolved in the water would have escaped from the 
water and would have to be vented. Above 205 °F the chamber would be 
pressurized to about 20 psig above the saturation pressure. Insulation 
was placed over the front viewing window to minimize the three-dimen¬ 
sional convection in the cavity. 

As the water in the chamber was heated to operating conditions, 
the clarity of the water was checked. If the chamber was not flushed 
adequately or if boiling occurred, there would occur a cloudiness which 
would become so dense that photographs could not be taken. In such a 
case the test would be abandoned. 

Temperature recordings were taken every half hour prior to nd 
throughout the photographic session. Photographs were taken when it 
appeared from the temperature recordings that the conditions were rela¬ 
tively stable, that is a drift of less than one degree per hour. 
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The camera was positioned in front of the viewing port and the 
trace particles were injected. The injector pressure had to be above 
the chamber pressure. The trace-injector valve was opened slightly to 
permit a small steady stream of particles to enter the modeling area. 
Great care was taken to regulate the rate of entry. If too much trace 
had been permitted to enter, the glass on top of the illuminated area 
might very well have been entirely covered with aluminum, blocking the 
light from above and forcing abandonment of the test. The trace was 
injected until it was sufficiently distributed throughout the chamber. 

The camera was focused on the particles in the illuminated area. 

To obtain the sharpest image the camera was first focused at widest 
aperture then stopped down to normal shooting aperture. Sufficient 
time (no less than 25 minutes) was allowed for the disturbance caused 
by the trace injection to die out before photographs were taken. Ex¬ 
posure time, light intensity, and aperture settings for photographs 
were determined by trial and error. Variation of these three parameters 
usually produced at least three out of ten usable exposures. The lights 
were raised to maximum intensity as quickly as possible during the last 
second of the exposure. This caused a bright dot to appear on one end 
of the trace image, and thus direction of flow could be denoted. The 
time was recorded at the completion of each exposure. 
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The film was developed as soon as possible after making the ex¬ 
posures in order to preserve the sharpest image. A high acutanee de¬ 
veloper (Rodinawas used because the softer working developers do not 
produce streaks of sufficient definition. A high-contrast grade of 
printing paper Agfa Brovia No. 6) was used for best results, and the 
paper developer tOaktol) was used in concentrated form to further in¬ 
crease the conLrast. 

About five hours were usually required to complete a run at one 
test condition. 


3. Experimental Results 


The basic convection producing mechanisms which influenced the 
heated water wore (1) r emperature difference between opposite walls 
of the model, ( 2 ) a simulated occupant whose temperature was above 
that of the chamber, and (3) natural ventilation due to a temperature 
difference between the water outside and inside the chamber. The tests 
have been segregated into four basic conditions which are as follows: 

Condition I : A chamber without ventilation or simulated 
occupant 

Condition II : A chamber with natural ventilation but 
without simulated occupant (s) 

Condition 111: A chamber without natural ventilation but 
with simulated occupant(s) 

Condition 1 ' : A chamber witii natural ventilation and 
with simulated occupant (s) 
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The results of the tests are presented in the following format: 
a narrative description of the fluid motion, a photograph of the flow, 
a schematic of the flow pattern shown in the photograph, and a figure 
showing the temperatures at specific locations recorded during the test. 

In all of the photographs of the flow patterns four vertical and 
four horizontal lines appear. This was the wire grid used to support 
the thermocouples. The grid was located about four inches behind the 
illuminated area; one thermocouple was fastened to each intersection. 

The thermocouples extended beyond the grid to a location immediately 
to the rear of the illuminated area. In many photographs two bright 
lines appear to be going into the ceiling. These were the thermocouples 
used to measure the ventilation supply and return temperatures. 

All temperatures shown on the following figures are in degrees 
Fahrenheit. 


Condition I 

Natural convection was induced only by a temperature difference 
between the walls, the right wall was maintained at a warmer temperature 
than the left, giving a counter-clockwise rotation in the modeling area. 
Teat results representative of Condition I are shown in figures 18-23. 
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The flow patterns (figures 19 and 22) due to the .XT between the 
walls were dependent on the water temperature and the magnitude of the 
.XT between the walls. Two water temperatures were considered for this 
situation--124 °F and 250 °F. A counter-clockwise vortex appeared in 
the bottom right corner at both of these temperatures. The flow from 
bottom to top occurred entirely in the right wall boundary layer. The 
interior flow pattern of the 250 °F test consisted mostly of horizontal 
motion while the interior flow pattern for the 124 °F test was primarily 
vertical from top to bottom. For the 250 ”F case this flow appeared to 
originate in the downward-moving left wall boundary layer and terminated 
near the right wall. A vortex existed in the upper left corner. A 
stream originating a half inch below the ceiling in the left wall 
boundary layer flowed toward the right wall. This stream turned down¬ 
ward just before touching the right wall boundary layer and terminated 
in the floor boundary layer. 


Condition II 

Ventilation ducts were introduced into the model to produce added 
convection effects. The simplest duct configuration consisted of ducts 
at the top right and top left of the chamber. The water flowed out of 
the chamber through one duct, was cooled, and then returned via the 
other duct. Test results representative of Condition II arc shown in 


figures 24-41. 
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Figure 22 Schematic of flow distribution for Condition 
I (refer to figure 21) 
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A characteristic flow pattern showing the streamlines for the sim¬ 
ple duct case can be seen in figure 25. Three-dimensional effects, i.e., 
water flowing out of the lighted two-dimensional slice into the surround¬ 
ing darkened water and vise-versa, were observed in these photographs. 
These three-dimensional effects appeared as sources and sinks in the 
photographs. In the flow-pattern drawings the sources are depicted as 
circles or ovals (without arrow points) and the sinks are indicated by 
X's. The term source implies flow into the lighted area from the sur¬ 
rounding water; a sink implies water flowing out of the lighted area. 

With the addition of ventilation ducts the flow pattern became 
considerably more complicated than in the non-vented case. In this 
cas there was a 2 °F AT between the walls and a 15 °F AT between the 
water entering the modeling area through one duct and leaving through 
the other duct. The water entering the model flowed down along the 
left wall to the floor. The water then flowed out across the floor, 
much of it passing out of the lighted field. There appeared to be a 
source near the right wall approximately a third of the distance from 
the top of the chamber, and a sink appeared in the bottom left corner. 
Most of the water that left the model through the exhaust duct did not 
come from the lighted field. The water entering the exhaust duct was 
supplied by a source located at the top center of the cavity. Some of 
the water entering the exhaust duct was supplied by the boundary-layer 
flow coming up the right wall and crossing the ceiling. 
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Schematic of flow distribution for Condition 
II (refer to figure 24) 
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Figure 26 Temperature distribution for Condition II 
(refer to figure 24) 
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Longer ventilation ducts were installed, extending approximately 
one inch into the cavity. A series of three tests were conducted with 
these new ducts. The temperature conditions for these tests are shown 
in the following table. 

Table No. 1 


Test 

No. 

Water Temp. 
(°F) 

between ventilation 
ducts 
(°F) 

AT between 
walls 
( °F) 

I 

250 

24 

4 

II 

250 

25 

2 

III 

300 

40 

2 


Figure 28 is the schematic flow diagram for Test I. The Clow was 
generally from right to left. There was a vortex in the bottom right 
comer and another one approximately two inches below the exhaust duct. 
Two sources existed: one immediately below and to the left of the ex¬ 
haust duct; the other, adjacent to the right side of the supply duct. 
There was a sink beneath the supply duct just above the floor. 
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Figure 27 Photo of convection occurring in chamber with 246 °F water, 
a 4 °F .XT between walls and a 22.6 XT between supply 
and exhaust ducts (Condition II, Test 1) 
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Schematic of flow distribution for Condition 
II, Test I (refer to figure 27) 
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Figure 29 Temperature distribution for Condition II, 
Test I (refer to figure 27) 
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The schematic flow diagram for Test II is shown in figure 31. The 
flow originates from the right wall boundary layer and a source beneath 
and to the left of the exhaust duct. The flow from the source did not 
move toward the exhaust duct as one might expect, but moved directly 
toward the left and joined the downward flowing water leaving the supply 
duct. Several vortices were present in this test. These were located 
in the bottom right comer, at the base of the water column leaving the 
supply duct, and at the ceiling midway between the two ducts. All of 
these vortices rotated in a counter-clockwise direction. 

Figure 34 is the flow diagram for Test III. The flow patterns in 
this test subdivided into three areas. In the upper third of the model 
the flow was generally from left to right. The streamlines in this 
upper third eventually ended at the exhaust duct. In the bottom left 
third of the model the flow was basically from right to left. These 
streams terminated in the downward flowing water from the supply duct. 
The flow in the bottom right third consisted of four vortices. These 
vortices were vertically stacked. The bottom vortex turned counter¬ 
clockwise. Each adjacent vortex had an opposite direction of spin. 

Another duct configuration simply extended the length of the 
return duct. This added length brought the bottom of this duct to 
within an inch of the floor. 
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Figure 30 ’hoto of convection in chamber with 250 "F water, a 2 I 
between wails and a 24 "F 1 between supply and exhaust 
(Condition II, Test II) 
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Schematic of flow distribution for Condition I 
Test II (refer to figure 30) 


































Figure 35 Temperature distribution for Condition II, 
Test III (refer to figure 33) 
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The flow shown in figure 37 was primarily from left to right. The 
stream originated in the left wall boundary layer and terminated in the 
right wall boundary layer. Deviations from this left to right flow oc¬ 
curred along the floor and ceiling. Approximately three quarters of an 
inch above the floor there was a stream moving from right to left. This 
stream was connected to the floor boundary layer by several downward 
moving streams. The upper third of the cavity mirrored the flow just 
described at the floor of the cavity. The volume of water moving in 
the upper third of the cavity was considerably greater than that movim: 
near the floor. 

One test was conducted with the ducts in the i'<«i fi gurat ion shown 
in figure 39. Right angle elbows, aimed toward the center of the cavity, 
were attached to the ends of the ducts. The average water temperature 
was 251 °F. The 'JC between the wall surfaces was 2.5 ’F, and between 
the convection ducts it was 15 °F. The flow pattern is shown in figure 
40. The streams move primarily from top to bottom. Approximately one 
inch above the floor a stream which moved horizontally from left to 
right picked up the vertical stream. This horizontal stream terminated 
in the left wall boundary layer. In the upper quarter of the cavity 
most of the streams originated in the ceiling boundary layer and, 
making a "U" turn, terminated in the exhaust duct. 








Figure 36 Photo of convection occurring in chamber with 196 °F water, a 
5 °F XT between walls and a 24.5 °F '\T between supply and 
exhaust ducts (Condition II) 
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Figure 38 Temperature distribution lor Condition II 
(refer to figure 36) 
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e 40 Schematic of flow distribution for Condition 
II (refer to figure 39) 
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Condition III 


This set of tests incorporated the use of one and then two simulated 
occupants. These occupants added heat to their surroundings and induced 
convection. Cylindrical figures were used specifying the size and heat 
dissipation of the simulated occupants. Test results representative of 
Condition III are shown in figures 42-47. 

The first of these tests was run with one occupant. The ventilation 
ducts were removed from the cavity. The left wall was kept 1.5 °F cooler 
than the right wall, and the average water temperature was 255 °F. The 
flow pattern for this particular case can be seen in figure 43. 

The -XT between the wall surfaces set up a counter-clockwise boundary 
layer flow. The simulated occupant, located at the center of the floor, 
caused a vertical upward flow of relatively high velocity. This stream 
divided to the right and left at the ceiling. The stream to the left of 
the occupant wandered in an "S"-like fashion toward the left wall boundary 
layer. A counter-clockwise turning vortex was formed just left of center 
at the ceiling. The flows off from the right of the occupant followed 
four different paths before reaching the boundary layer. The uppermost 
started a left to right boundary layer flow. This flow was in opposition 
to the existing boundary layer flow and resulted in turbulence along 
the ceiling. The second highest flow turned downward and to the right. 
This intersected the right wall boundary layer approximately two inches 
from the ceiling. The third highest flow followed an "S" course to the 
floor boundary layer. The lowest flow followed a similar "S" path to 
within an inch of the floor. At this point the flow turned to the left 
and united with the floor boundary layer approximately one incti from the 
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Schematic of flow distribution for Condit 
III (refer to figure 42) 





















left wall 


One observation of this single occupant test is particularly note¬ 
worthy. Except boundary layer flows, all flows that move from top to 
bottom move toward the center of the cavity as they pass by the top of 
the occupant. After passing the top of the occupant, they move outward 
again toward the sides of the cavity. 

For the next test involving two simulated occupants, 250 °F water 
was used and there was a 2.5 °F between the inside wall surfaces. 

No convection ducts were used for this test. The resulting flow pattern 
cai- I'c seen in figure 46. 

Three-dimensional effects were prevalent in this two-occupant test. 
Streamlines commenced and terminated without appearing to have a logical 
source or sink. A vortex occurred in the upper left corner. This vor¬ 
tex was fed primarily by the ceiling boundary layer. A flow leaving 
the left wall boundary layer, just below this vortex, traveled downward 
and to the right. This flow bifurcated just prior to touching the verti¬ 
cal flow leaving the left occupant. One branch mated with the flow leaving 
the occupant. The other branch flowed downward and to the left, termi¬ 
nating in the left wall boundary layer. Another flow appeared to origi¬ 
nate from behind the left occupant about a half inch from the top of the 
occupant. This flow moved to the left, turned downward, and then flowed 
back to the right whereupon it divided into two branches before disap¬ 
pearing behind the left occupant. 
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Figure 45 Photo of convection occurring in chamber with 242 °F water 
and a 3 °F XT between walls (two simulated occupants. Con¬ 
dition III) 
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Schematic of flow distribution for Condition 
III (refer to figure 45) 











The flow along the ceiling was turbulent. Half way between the 
occupants, approximately one inch from the coiling, a source fed a flow 
that terminated in the vertically rising currents of both occupants. 

There was considerable turbulence and vorticity along the ceiling to 
the right of the right occupant. Two skewed concentric flow loops 
originated and terminated on the right side of the rising currents 
leaving the right occupant. In the bottom right corner was a vortex 
which turned counter-clockwise. A flow originated in the right wall 
boundary layer and terminated behind the right occupant. 

Condition IV 

Natural convection was induced by a combination of the temperature 
difference between the walls, the presence of simulated occupants, and 
natural ventilation to the modeling area. Test results representative 
of Condition IV can be seen in figures 48-53. 

Ventilation ducts were used in this test which involved one simu¬ 
lated occupant. These ducts did not penetrate the ceiling of the cavity, 
i.e., the openings were flush with the ceiling. The water temperature 
was 250 °F, the between the wall surfaces was 1 °F, and the XT be¬ 
tween the convection ducts was 30 °F. The flow pattern for this test 
can be seen in figure 49. 
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Figure 48 Photo of convection occurring in chamber with 247 °F water, a 
2,1 °F I between wails and a 27 °F \T between supply and ex¬ 
haust ducts (one simulated occupant, Condition IV) 
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Figure 50 1emporaturt distribution lor Condit 

(refer to figure 48) 
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The flow in this case, as in the previous one-occupant case, was 
divided down the middle by the rising currents leaving the occupant. 

Upon reaching the ceiling the flow from the left side of the occupant 
moved toward the left. The flow moved downward and to the left until 
it united with the columnated downward flowing water leaving the supply 
duct. Two counter-clockwise turning vortexes were located just to the 
top left and bottom left of the occupant. The flow leaving the right 
side of the occupant rose to the ceiling. After arriving at the ceil¬ 
ing the stream moved to the right toward the exhaust duct. Consider¬ 
able turbulence occurred at the ceiling. A clock-wise-rotating vortex 
was formed approximately one and a half inches below the exhaust duct. 
This vortex was primarily fed by a stream leaving the vertical water 
column above the occupant. A second vortex, turning counter-clockwise, 
was formed about one inch beneath the first vortex. Downward-flowing 
streams connected the two vortexes. A stream leaving the boundary 
layer near the bottom right corner bent back and fed both the bottom 
vortex and the upward-flowing boundary layer on the right side of the 
occupant. 

The second test for this condition involved two occupants. The 
water temperature was 250 °F, the .XT between the wall surfaces was 3 °F, 
and the XT between the two ducts was 12 °F. The flow pattern for this 
test can be seen in figure 52. 
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In this test, as in the previous two-occupant test, three-dimen¬ 
sional effe cs were observed. The velocity of the water entering and 
leaving the ducts was considerably slower than in any of the previous 
tests. This low velocity was due to the relatively small maintained 
between the two ducts. 

The rising currents from the right occupant fed both the ceiling 
boundary layer and the exhaust duct. The ceiling boundary layer moved 
from right to left. Two streams, originating from the rising currents 
of the right occupant, flowed horizontally to the right wall boundary 
layer. The lower of these two streams divided into two branches and one 
of these branches terminated in the boundary layer. The second of these 
branches flowed behind the occupant. A vortex, turning counter-clockwise, 
was established in the bottom right corner. All the flow between the 
occupants appeared to originate from a source behind the rising currents 
of the left occupant. One stream from this source flowed upward and 
ended in the ceiling boundary layer. Another flowed downward and termi¬ 
nated in the floor boundary layer. Both of these streams bifurcated. 

The two resulting secondary streams flowed behind the right occupant. 

The water leaving the supply duct flowed in an "S" path to the floor 
boundary layer. One stream from the supply duct and another from the 
rising currents of the left occupant combined to produce a vortex which 
rotated counter-clockwise. This vortex was approximately two inches 


above the left occupant. 









6. Discussion of Method and Results 


The tracer particles used in this study were very small flakes oi 
aluminum. Since the density of aluminum is greater than that of water, 
these flakes eventual';, settled at the bottom of tile chamber. At first 
it was expected mat some discrepancy would exist between the motion 
indicated by the traces on the photograph and the true flu-d mot ion. 

The settling velocity, however, was found to be quite small. For ex¬ 
ample, when the tracer was put in room temperature water (72 °F), most 
of it remained suspended in the water overnight. When put in warm 
water (120 °F) most ol the trace remained suspended for four hours or 
more, and when put in JOC T water about 50. of the original amount 
settled in about one half tour. 

Due to the shape cl the aluminum particle there was room for error 
in the velocity measurements. The particle is flat in one plane and a 
multi-sided polygon when viewed in a plane perpendicular to the first. 
When the flake moved in 11 to fluid it tended to align with the- direction 
of shear of the fluid. If tile flow were laminar, a flake could move a 
considerable distance exposing only one side eo the viewer. If the 
fluid motion were turbulent, the particle’ tended to tumble. A particle 
coul1 present an edge view to the film and leave no trace at all. When 
the traces made by the flakes wore used to determine the velocity of the 
fluid it was often difficult to tell whether the length of the streak 
indicated the actual distance movi-d or just the distance moved when 
the trace was reflecting light to the iilm. Spherical particles would 
probably be the most desirable shape for velocity measurement hi cause 
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they would constantly reflect the same amount of light regardless of 
orientation. But of those tried none were of sufficient density to re¬ 
main suspended in the high-temperature water. To reduce the chance of 
error due to the particle shape, a fairly large number of streaks were 
measured for each velocity determined. 

The model was essentially designed for a study of two-dimensional 
convection. The apparatus was constructed with sufficient depth to 
minimize end effects at the midsection of the model enclosure. While 
the end effects were reduced, there were still three-dimensional effects 
due to a longitudinal temperature gradient. When the glass plate viewing 
port was insulated on its outside surface, it was usually two or three 
degrees cooler than the rear wall and the three-dimensional effects were 
small. But when the insulation was removed for fifteen or twenty minutes 
during photographic sessions, the temperature difference between front 
and rear became greater and these three-dimensional effects increased. 

The error due to this factor was minimized by taking the photographs 
immediately after the insulation was removed from the glass viewing 
port. 
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The Prandtl number of room temperature air is 0.72. Water at 
500 °F has approximately the same Prandtl number as room air (0.87). 
Further, a Grashof number equivalent to that found in a full-sized 
room can be obtained by using high temperature water in a small-scale 
model. Another attractive feature of high temperature water is that 
above 350 °F the Prandtl number is relatively constant with respect to 
temperature change. This characteristic minimizes inaccuracies due to 
property value changes across the thermal boundary layers next to the 
heat transfer surfaces. 

Because of physical limitations of the modeling apparatus used 
here, water at 250 °F was used for the simulations. The saturation 
pressure of water at 250 °F is about 30 psia, while at 500 °F it is 
about 680 psia, thus a considerably stronger container is required 
for higher temperature simulations. Also rusting occurred in this 
particular apparatus at temperatures above 250 °F. As the temperature 
increased rusting occurred to such an extent that the trace particles 
in the water were obscured. 

By using 251 °F for this series of simulations, the same Grashof 
number, and aspect ratio, and ncarl the same Prandtl number (Gr — 10 
L = 1., Pr = 1.A1) were reproduced in the small-scale model (6" x b ,! , 
width x height) as are normally found (Gr — 10^, L = 1., Pr = 0.72) 
in a full sized room (8' ;• S', width x height). 
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One test was run at an average water temperature of 124 °F to show 
the difference between convection at this temperature and the other tests 
of 250 °F. Figures 18, 19, and 20 show the results of that test. In 
the low temperature test the flow is slower, more viscous, and there is 
much less turbulence. The Grashof number is an order of magnitude less 
for the low temperature test. 

The differences between the two tests are due almost entirely to 
the viscosity of the water., The viscosity of the 124° water is nearly 
twice that of the 250° water. The lower the viscosity the higher the 
Grashof number and the taster the flow duo to a given temperature dif¬ 
ference . 

Figures 33, 34, and 35 show the results of a simulation at an 
average water temperature of 292 °F, approximately 40 degrees above 
normal simulation temperatures. At this higher temperature the Prandtl 
number of water (about 1.1) more closely approximates the Prandtl number 
of air (0.72). Thus operation at this level is more desirable for an 
accurate simulation. However, as noted, there were considerable ex¬ 
perimental difficulties involved with operating at this high temperature. 
The fact that results were obtained at all was largely due to good 
fortune. 
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The effect of the increased temperature on predicted fluid motion 
is apparent. Since the effect of the higher temperature is to lower tin 


viscosity of the water, the same temperature differentials across the 
model operated at a higher temperature level would have high velocities. 
The higher velocities would increase the ventilation rate and turbulence 
level. Further, when the ventilation rate is high the temperature dif¬ 
ference between the floor and ceiling is smaller as was observed here. 

A comparison of the results obtained from small-scale modeling of 
test Conditions I - IV with the results of similar tests in a full-sized 
model would certainly be desirable. Unfortunately there has been little 
previous work in full-scale models that can be compared quantitativelv 
with the results of these simulations. A comparison cannot be made 
because the boundary conditions of other full scale tests were not the 
same as those in this model, or because the boundary conditions were 
not available in the reports of the other tests. 

However, a qualitative comparison with Daiv/s ' results [15] may be 
made with some interesting agreements. Daws observed that for any flow 
of air moving along a wall a point downstream is reached where the 
velocities of tii it flow are comparable with those of the room. At 
tiiat point, which is dependent upon the strengtli of the flow, all en¬ 
trainment ceases and the stream spreads out in all directions, breaking 
away from the wall and joining movement induced by entrainment elsewhere'. 
This phenomena llso was observed in this small-scale modeling study. 

When the temperature differential between the ventilation supply and 
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The effect of the increased temperature on predicted fluid motion 
is apparent. Since the effect of the higher temperature is to lower tin 
viscosity of the water, the same temperature differentials across the 
model operated at a higher temperature level wouid have high velocities. 
The higher velocities would increase the ventilation rate and turbulence 
level. Further, when the ventilation rate is high the temperature dif¬ 
ference between the fLoor and ceiling is smaller as was observed her'. 

A comparison of the results obtained from small-scale modeling of 
test Conditions I - IV with the results of similar tests in a full-sized 
model would certainly be desirable. Unfortunately there lias been little 
previous work in full-scale models that can be compared quantitatively 
with the results of these simulations. A comparison cannot be made 
because the boundary conditions of other full scale tests were not the 
same as those in this model, or because the boundary conditions were 
not available in the reports of the other tests. 

However, a qualitative comparison with Daws' results [15] may be 
made with some interesting agreements. Daws observed that for any flow 
of air moving along a wall a point downstream is reached where the 
velocities of ti> u flow are comparable ’.ith those of the room. At 
that point, which is dependent upon the strength of the flow, all en¬ 
trainment ceases and the stream spreads out in all directions, breaking 
away from the wall and joining movement induced by entrainment elsewho-e. 
Tiiis phenomena ilso was observed in this smalt-sca modeling study. 


When the temperature differential between the ventilation supply and 






exhaust air was small the supply air was observed to move down the left 
wall for a short distance, then to mix with the room as described above. 

Daws noted that air streams close to a room surface are drawn to 
the surface by what appears to be a slight sucking action at the wall 
surface (known as the Coanda effect). This phenomena was also observed 
in the small-scale model when the ventilation streams entered the model. 
As the ventilation supply air moved toward the floor it was drawn to 
the wall. The distance from the ceiling to the point where the fresh 
air stream came into contact with the wall was dependent upon the 
strength of the stream. The stronger the flow the more it resisted 
the sucking effect of the wall. 

Daws noted that the action of the air leaving the room was confined 
to the immediate vicinity of the outlet. This was also observed in the 
small-scale model simulations. 

The photographs in Daws’ paper show the horizontal motion of room 
air due to entrainment with vertically moving streams. This occurrance 
can also be seen in the photographs of small-scale model simulations. 
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7. Summary 


The problem o£ obtaining a quantitative velocity distribution map 
of a natural convection field has become a classic of engineering re¬ 
search. Numerous approaches to the problem have been made using botli 
theoretical and experimental techniques. While the experimental method 
presented here is not the complete answer, it is a step toward solution 
of the problem. Limitations within the method as well as with this 
particular experimental apparatus prevent it from becoming the complete 
answer. 

The experimental apparatus used here is suitable for two-dimensional 
flow studies only. An apparatus can be designed and constructed to ex¬ 
tend the method to three-dimensional studies. 

The experiment discussed in this report is not conclusive proof 
that the method is qualitatively accurate, even for a two-dimensional 
simulation. A comparison with the results of a dynamically similar 
full scale test is required. For example, it is necessary to know the 
difference between these experimental results and full scale results 
as a function of the Prandtl number. It may be that the overall air 
ilow pattern, which is the dominating influence on the occupants' com¬ 
fort, is not as dependent on the Prandtl number as the boundary layer 
velocity profile Another factor leading to inaccuracies between the 
results observed in this model and those in a prototype room is the 
dissimilarity of boundary conditions. The temperature distribution 
over the walls of this model do not closely represent those over the 
surface of an actual room wall. Greater control of wall temperature 
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distribution must be built into a future model. The magnitude of the^e 
effects are not known and it can only be determined by tests in a full 
sized room. 

In spite of these shortcomings in rigorous duplication of actual 
room conditions, it is strongly felt that an ultimate solution to 
natural convection velocity patterns may be attained through modeling 
techniques. The construction advantages of a model over a full scale 
apparatus are more than enough to merit further development of this 


technique. 
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Appendix 


Simulated Air Distribution in a Single Room Shelter 

The results presented in the main body ol Lhis report have been 
transformed into a simulated case of air in a full sized shelter and 
are presented in this appendix. Quantified data, such as temperatures 
and velocities, have been scaled through the appropriate similarity 
parameters. Flow patterns shown here are identical to those presented 
in Chapter 5. The order of presentation and the types of conditions 
considered also follow those in Chapter 5. 





Condition I: Simulation of Air Distribution in a 
Shelter Without Ventilation or Simulated Occupant 

This condition represents the simplest set of boundary conditions 
considered. Figures A1 and A2 show the air current patterns and tempera¬ 
ture and velocity profiles predicted by the small-scale model for the 
full-scale shelter. The size of the full-scale shelter is 8 feet wide 
by 8 feet high, the length is not considered. 

Several features are apparent from observing the predicted air tem¬ 
peratures. The first is the thermal stratification of the air at the 
inner core of the room. The inner core being defined as that volume 
of air extending from the center of the room to within about one foot 
of all the surfaces. This stratification may be observed throughout 
the simulations for all four conditions studied. 

The cause for the flow shown in figure Al, is the preset temperature 
difference between the left and right walls. Since the right wall is 
warmer than the left a thin boundary layer will flow up the right wall, 
cross the ceiling and then fall down the cooler left wall to the floor. 

On reaching the floor the flow will continue to hug the boundary while 
crossing to the right wall. This produces a general counter-clockwise 
rotation at the interface of the air and the room surfaces. It also 
produces a secondary motion in the inner core air. The secondary motion 
is due to the inner air becoming entrained with the rapid moving boundary 
layers along the walls. 
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Air distribution predicted tor an 8' x 8' shelter 
without natural ventilation simulated occupant 
(Condition I) 
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A quantitative estimate of tlie boundary layer and its velocity in 
tlie model is not noted due to experimental limitations. It is suggested 
though, to be several times thinner and faster than the flows across the 
ceiling and floors. Frcm the measurements made, the predicted air flow 
across the ceiling is approximately 3.5 feet-per-minute, and the boundary 
layer would be on the order of 9 in. to 12 in. thick. The predicted air 
flow across the floor, while in the opposite direction is about tlie same 
velocity and thickness. The predicted velocity of the air in the inner 
core is considerably lower and mainly in a horizontal direction. The 
maximum core air velocity would be near the center of the room and 
slightly over one foot per minute. 

Condition II: Simulation of Air Distribution in a Shelter 
With Ventilation, But Without a Simulated Occupant 

Simulations were made for a number of variations on this condition 
The boundary conditions varied were the temperature difference between 
the left and right walls, the temperature difference between the ventila¬ 
tion supply and exhaust air, the temperature level of the modeling fluid 
and the positions of the ventilation duct openings. 

Figures A3 and A4 show the air distribution patterns and temperature 
and velocity profiles predicted for a full sized (8 1 x 8') naturally 
ventilated shelter. The ventilation ducts are flush with the ceiling 
and in a full sized shelter would be located approximately 18 in. from 
either side wan. The inside diameter of the ventilation ducts are 
equivalent in size to 3.75 in. round ducts. 
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Figure A3 Air distribution predicted for an 8' x 8' 
shelter with natural ventilation but 
without a simulated occupant (Condition II) 
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The same thermal stratification of the inner core observed in Con¬ 
dition I can be seen here. Assume a one degree temperature difference 
exists between the left and right walls, while a comparatively large 
difference (15 °F) exists between the ventilation supply and exhaust air. 
this indicates that most of the motion is due to the ventilation. The 
predicted temperature of the air leaving through the ventilation exhaust 
is slightly above the predicted temperature of the air near the ceiling. 
This is due to some of the warmer air from the right wail boundary layer 
mixing with the warm air near the ceiling as both were exhausted. 

The predicted velocities of the different air flows in the full-sized 
shelter are shown in figure A4. The cold incoming air falls directly 
to the bottom of the shelter. On its way to the floor it is drawn towards 
the left wall, due to the Coanda effect. The entering air entrains the 
room air as it falls, and causes a horizontal motion of the room air 
toward the incoming stream. The horizontal velocity of the room air is 
highest adjacent to the opening to the ventilation supply duct, since 
the entrainment of room air with the entering air is highest at this 
location. 

As the ventilation air moves down the wall its velocity diminishes 
and the fresh air becomes increasingly mixed with the room air. Upon 
reaching the floor the ventilation air spreads out across the floor. 
Although a little of the fresh air may become entrained in the slight 
boundary layer rising along the right wall, most of it remains in the 
lower portion of the room. 
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The velocity of the incoming ventilation air has been determined 
approximately three feet below the supply duct opening. The maximum 
predicted velocity of the air stream for this case is 29.5 feet-per- 
minute. The horizontal velocity near the center of the room would be 
1.1 feet-per-minute and the maximum velocity of the flow across the 
floor would be 10.2 fect-per-minute. The boundary layer at the floor 
would be about 9 in. thick. 

Several variations on the above simulation have been made for dif¬ 
ferent temperature differences between the ventilation supply and exhaust 
ducts. Results of these simulations will be discussed, but the figures 
will not be shown since the shelter configurations are the same as those 
already discussed. 

If the temperature differential between the ducts is increased to 
29 degrees, considerably more turbulence will occur in the shelter. The 
room air will mix with the fresh air much more rapidly, the horizontal 
velocities due to entrainment with the supply air will be higher than 
before, and the temperature differential between the floor and ceiling 
will decrease. The velocity of the entering air, measured at the same 
location as before, will have a maximum value of 41 feet-per-minute. 

The temperature difference between the floor and ceiling will be reduced 
to 4 °F. 

If the temperature difference between the supply and exhaust is re¬ 
duced to six degrees, a significant reduction in the air motion will 
occur throughout the room. As the air enters it will be deflected toward 
the left wall close to the inlet. The fresh air will move down the wall 


a short distance and then turn horizontal about halfway down the wall. 








The velocity of the supply air, at the same location as befor-', will be 
a maximum of 17.8 feet-per-minute. The temperature differential between 
the floor and ceiling will be 9° higher than for the original case shown 
in figure A4. 

The next set of figures depicts results for some other tests conducted 
under Condition II. The significant difference between these tests and 
those already discussed for Condition II is the position of the ventilation 
ducts. 

Figures A5 and A6 show the combined effect on the room air, of a 23 
degree F temperature difference between the supply and exhaust air and 
a 5 degree difference between the left and right wall temperatures. 

Here the ventilation ducts are protruding down about three feet from the 
ceiling. The considerable horizontal motion of the room air is due to 
its entrainment with the faster, vertically moving streams. Since the 
supply air enters the room closer to the floor it undergoes less mixing 
with the room air before reaching the floor. After spreading out across 
the floor there is little further mixing with the room air, except for 
the small amount that is entrained in the boundary layer rising along 
the right wall. Due to the low position of the opening to the exhaust 
duct, air is removed from that elevation. This causes a layer of warm 
stale air to build up near the ceiling. 

Figures A7 and A8 show the results of a simulation closely resembling 
the one just discussed. The significant difference here is the small, 
two degree, temperature difference between the left and right walls. The 
ventilation rate is almost exactly the same as before. Thus, the reduced 
air motion throughout the room is due to the small circulation effect of 


97 






tlie boundary layer. The net effect of the smaller temperature difference 
between the walls is to cause a larger temperature difference between the 
floor and ceiling. 

Ir. figures A9 and A10 can be seen the effects of placing the ventila¬ 
tion supply duct close to the floor. In a full-sized room this would be 
equivalent to the duct opening about 18 ir above the floor. Mixing is 
poor because the air enters close, to the floor, then spreads out in all 
directions at the floor level. Some circulation occurs due to the 5° 
differential between the left and right walls, but this circulation does 
not cause significant mixing of the air in the shelter. The relatively 
large, 10 degree, temperature difference between the floor and ceiling 
is further evidence of the poor mixing. 

Although the mixing is poor, this simulation indicates a definite 
advantage to be gained from locating the ventilation supply close to the 
floor of a shelter, rather than close to the ceiling. When the temperature 
differential between the supply and exhaust air becomes small, as it may 
during the summer months, the ventilation rate will be small. As has 
already been shown, when the ventilation rate is small and the supply 
duct opening is close to the ceiling, the entering air may never reach 
the shelter occupants. On the other hand, if the ventilation supply duct 
opening is located near the floor then the occupants will benefit from 
even a small amount of fresh air. 

In figures All and A12 can be seen the results of a simulation when 
the direction of the openings to the ventilation ducts is changed from 
vertical to horizontal. One of the first things to notice is the apparent 
good mixing of the room air with the ventilation air, in spite of the 
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small, 14 degree, difference between supply and exhaust temperatures. 

This mixing is further evidenced by the relatively small, 5 degree, tem¬ 
perature difference between the floor and ceiling. 

Actually the mixing throughout the shelLer is not quite as good as 
it appears when one considers a three-dimensional shelter. The supply 
duct is highly directional thus the ventilation air is discharged directly 
toward the right wall. Once the flow reaches the right wall it does not 
spread out uniformly across the floor but tends to divide and go toward 
the rear and front of the shelter. The temperature difference between 
the walls was relatively small; thus, it did not contribute in any great 
part to the circulation in the shelter. 

The advantage of having the ventilation duct discharge horizontally 
is that the flow of air can be directed toward a specific location in 
the shelter. This would be a disadvantage however, in a crowded shelter, 
because many of the shelter occupants would be missed by this highly 
directional stream. 

The changed direction of the exhaust duct opening had no apparent 
effect upon the air distribution in the shelter. 

Condition III: Simulation of Air Distribution in a Shelter 
Without Natural Ventilation But With a Simulated Occupant 

Condition III resembles the situation that would occur when the 
outside air temperature exceeded the shelter temperature or when the 
shelter must be operated under sealed conditions. 

Figures A13 and A14 show predicted air current patterns and velocitie 
and temperatures in a full-sized, S fer- - J de by 8 feet high, naturally 


ventilated shelter 
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e A5 Air distribution predicted for an 8' x 8' shelter 
with natural ventilation but without a simulated 
occupant (Condition II) 


100 








64 


Figure A6 Air temperatures predicted for an 8' x 8' 

shelter with natural ventilation but without 
a simulated occupant (Condition II) 
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Air distribution predicted for an 8' x 8' 
shelter with natural ventilation but without 
a simulated occupant (Condition II) 
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Figure A8 Air temperatures predicted for an 8' x 8' shelter 
with natural ventilation but without a simulated 
occupant (Condition II) 
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Air distribution predicted for an 8 x 8* 
shelter with natural ventilation but without 
a simulated occupant (Condition II) 
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Figure AlO Air temperatures predicted for an 8' x 8' 
shelter, with natural ventilation but 
without a simulated occupant (Condition II) 
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All Air distribution predicted for an 8' x 8' 

shelter with natural ventilation but without 
a simulated occupant (Condition XI) 
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Figure A13 Air distribution predicted for an 8' x 8' 

shelter with a simulated occupant but without 
natural ventilation (Condition III) 
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Figure A14 Air temperatures and velocities predicted 
for an 8' x 8' shelter with a simulated 
occupant (77 °F) but without natural 
ventilation (Condition III) 
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The temperature of the occupant*s surface is about 17 °F above the 
floor temperature. The velocity of the air rising from the occupant would 
be 31.2 feet-per-minute. 

One of the most apparent effects of Condition III is the large tem¬ 
perature differential, 22 degrees, between the floor and ceiling. In 
the absence of ventilation air to mix with and cool the room air, the 
heat from the simulated occupant is trapped at the ceiling. Further, a 
small, 1 degree, temperature difference between the left and right 
walls did not give a significant boost to circulation. 

In terms of an actual shelter it seems surprising that there should 
be such a large temperature differential between the floor and ceiling. 

It seems especially surprising that there would be a 17 degree temperature 
difference between the floor and the room air two feet above the floor. 

In general there is little motion in the shelter, except for that 
above the simulated occupant. The air rising from the occupant, despite 
its comparatively high velocity, does not produce a significant effect 
on the shelter air motion. The air warmed by the occupant moves vertically 
through the stagnant inner core until reaching the ceiling where it 
spreads out in all directions. As it passes through the inner core the 
rising column causes horizontal motion in the surrounding air. This 
happens because the surrounding air becomes entrained with the fast 
rising column of warm air. 

Figures A15 and A16 show the same shelter but with two simulated 
occupants. Again notice even a larger temperature difference between 
floor and ceiling. This is due to the presence of an additional simu¬ 
lated occupant. 
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Figure A15 Air distribution predicted for an 8' x 8' 

shelter with simulated occupants but without 
natural ventilation (Condition III) 


111 





















The highest velocities in the shelter occur above the simulated 
occupants. The warm air rises to the ceiling then spreads out in all 
directions. The temperature difference between the right and left walls 
is 3 degrees which is enough to cause some circulation around the periphery 
of the room, but not enough to have a significant effect on the mixing 
of the room air. 

Condition IV: Simulation of Air Distribution in a 
Shelter With Natural Ventilation and With a Simulated Occupant 

This is the condition most likely to be encountered in a shelter, 
provided there is a favorable thermal driving force between the inside 
and the outside. 

Figures A17 and A18 show predicted air current patterns and tem¬ 
peratures and velocity profiles for this condition. 

The simulated occupant is at a temperature approximately 17 degrees 
above that of the floor. The velocity of the air rising from the occupant 
is slightly over 30 feet-per-minute. 

Here there is roughly a 10 degree difference between the floor and 
ceiling whereas that difference is 22 degrees for Condition III. The 
ventilation effect would be primarily responsible for the smaller difference 
since it prevents the buildup near the ceiling of heat generated by the 
occupant. The temperature difference between the walls, 3 degrees, would 
cause a counter-clockwise circulation at the periphery of the air space, 
but it would not be the predominant influence on the air motion. 






A17 Air distribution predicted for an 8' x 
shelter with a simulated occupant and 
natural ventilation (Condition IV) 
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Figure A18 Air temperatures and velocities predicted 
for an 8' x 8' shelter with a simulated 
occupant and natural ventilation (Condi¬ 
tion IV) 
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There would be considerable air motion due to the combined effect 
of the heat from the occupant, the ventilation, and the circulation due 
to the temperature difference between : ,«_■ walis. 

Some thermal stratification noted in earlier cases would also be 
present here in spite of the rather vigorous air motion. This indicates 
that the inner core of air would be quite stable, and fairly resistant 
to the pert-orbing effect of streams passing through. Because of its 
large inertia in relation to that of the stream passing through, the 
core air would toad to absorb these streams or dissipate their effect, 
rather than to be excited by them. 

Figures A19 and A20 show a ventilated shelter with two simulated 
occupants. Again the temperature of both occupants is assumed to be 
about 17 degrees above the floor temperature. The velocity of the air 
rising from the occupants would be approximately the same as before, 
slightly over 30 feet-per-minute. 

The temperature difference between the floor and ceiling, about 
ten degrees, is somewhat higher than in the preceding simulation but 
this is due to the reduced ventilation rate. The heat rising from the 
occupants is trapped near the ceiling and drawn off gradually through 
the ventilation exhaust. Again some thermal stratification is apparent 
in this case. 
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A19 Air distribution predicted for an 8* x 8' 

shelter with simulated occupants and natural 
ventilation (Condition IV) 
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The ventilation effect would be small because the difference be¬ 
tween the supply and exhaust temperature is only 7 degrees. The air 
would enter the shelter and then be quickly drawn over toward the near 
wall. It would fall down the wall a short distance and then move hori¬ 
zontally toward the warm stream rising from the left occupant. The 
fresh air becomes entrained in the warm air rising from the left occupant 
and is returned to the ceiling. In this manner the fresh air would mix 
with the air near the upper part of the room and be eventually drawn 
off by the exhaust. This indicates, in effect, a "short circuiting" of 
ventilation air when the thermal driving force is small. 

Based on the progression of simulations presented, the separate 
effects of the different convection producing mechanisms can be under¬ 
stood. Ventilation has by far been the greatest influence on the air 
distribution. The air rising from a single occupant while of surprisingly 
high velocity does not have as significant an effect on the distribution 
of air in the room as does the ventilation air, principally because the 
air warmed by the occupant is not generally as large a volume as the 
ventilation air. In the case where there are a number of occupants in 
the shelter their combined effect would be significant. The least effect 
on air distribution would be that of the circulation due to the temperature 
difference between the walls. If the temperature differences between the 
walls is of the order of twenty to thirty degrees, then the effect might 
be more significant. However, the condition of a large temperature dif¬ 
ference between walls is not typical of a below ground shelter. 
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Conelusions 


A number of conclusions regarding air distribution in a naturally 
ventilated shelter may be drawn from the results of these approximate 
simulations. Some of the more significant conclusions follow: 

a. The temperature differential between the air inside and 
outside the shelter determines the rate of natural venti¬ 
lation to the shelter. 

b. A relatively large volume of the air in a naturally 
ventilated shelter is slow moving, has a surprising 
degree of thermal stratification, and is resistant to 
perturbing effects such as air currents produced by 
occupants or ventilation. This volume of air, referred 
to here as the inner core of air, extends outward from 
the center of the shelter to within about one foot of 
the walls, ceiling and floor. 

c. The inner core acts as a sink for the energy of turbulent 
flows. That energy is dissipated by vortices produced 

at the interface between the inner core air and fast 
moving flows which pass through the inner core. 

d. Vertical flows, within a foot or so of the wall, are 
drawn to the wall. This is due to the Coanda effect 
which gives the appearance of the wall drawing the 
flow toward it. 
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e. In those instances where a temperature difference 
exists between two opposite walls in a room or shelter, 
a flow of air is produced which circulates around the 
periphery of the room air space. The air rises along 
the warmer wall, crosses the ceiling to the cooler wall, 
drops along the cooler wall to the floor, and moves back 
across the floor to the warmer wall. The strength of 
the circulation is dependent upon the temperature dif¬ 
ferential between the two opposite walls. 

f. At the warmer wall a boundary layer occurs which flows 
quickly up the wall. The thickness of the boundary 
layer predicted for air is between one and two inches. 

At the ceiling the speed of the boundary layer decreases 
and its thickness increases. The width of the predicted 
air boundary layer across the ceiling is predicted to 

be about nine inches. As the air falls along the cooler 
wall it returns to a layer about one to two inches thick. 
The flow across the floor although opposite in direction 
is much like that across the ceiling both in thickness 
and speed. 
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g. The temperature difference between the floor and 
ceiling of the shelter is dependent on the extent 
of air mixing within the shifter. If the mixing 
is good, as caused by a high rate of ventilation, 
then the temperature difference will be small. 

When the mixing is poor, as when the ventilation 
rate is low then the air in the shelter tends to 
stratify and produce a rather large temperature dif¬ 
ferential between floor and ceiling. 

h. When the openings to the ventilation supply and return 
ducts are placed at the same elevation in a shelter 
and a low temperature differential is imposed between 
inside and outside air, then a short circuiting of 
the fresh air from the supply to the exhaust will 
occur. 

i. The optimum location for a natural ventilation supply 
duct outlet is close to the floor. This location 
minimizes the possibility of short circuiting between 
supply and exhaust, and also puts the fresh air close 
to the occupants. 

j. The optimum location for a natural ventilation exhaust 
duct outlet is at the ceiling level. The warmest stale 
air in the room is close to the ceiling and can be 
drawn o£i by an exhaust at that location. 
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k. The direction in which the opening to the exhaust duct 
faces had no influence on where the air is drawn from. 
Generally the air is drawn from a spherical area around 
the exhaust opening. 

l. The direction in which the supply duct discharges has a 
significant effect on where the air goes. 

m. When the temperature of the outside air equalled or ex¬ 
ceeded that of the air inside the shelter then all 
natural ventilation stopped. 

n. The maximum predicted velocity of air rising from an 
occupant was measured to be about 31 feet-per-minute 
and the maximum predicted velocity of the ventilation 
supply air, measured at a location three feet below 
the duct opening was about 41 feet-per-minute. Because 
the Prandtl number of the fluid was not closer to that 
of air, the velocities measured in the model are below 
those expected for air in a full-sized shelter. 
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Recommendations 


As a result of tiiese simulations a number of recommendations can bo 
made for the design of naturally ventilated below ground survival shel~ 
ters. These recommendations are as follows: 

a. All shelters should have manually operated ventilating 
fans to change the air within the shelter. During iiot 
days there is no assurance that the outside air tem¬ 
perature will not approach or even exceed the inside 
temperature. Once the temperature outside is higher 
than inside all natural ventilation ceases. 

b. The opening to the supply duct should be close to the 
floor rather than close to the ceiling. This discharges 
air into the shelter close to where the occupants are 
located and prevents short circuiting of the air between 
the supply and exhaust ducts. 

c. The supply air should discharge directly toward the floor 
rather than be directed out across the floor by a register. 

This causes the air to spread out across the floor more 
uniformly. 

d. A diffuser or register should not be put over the opening 
of the ventilation supply duct since it adds resistance 
to the flow of air. 

e. If the ventilation supply duct is inside the shelter, 

it should be insulated to minimise heating of the entering 
air. 
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f. The ventilation supply duct for new shelter designs 
could well be run through the earth outside the shelter. 
In this way the added cooling effect of the earth in 
the summer could be used to improve the thermal dif¬ 
ferential between the inside and outside of the shelter 
and hence to improve the rate of natural ventilation. 

In the winter the earth would be warmer than the ou side 
air and would tend to warm it to a more comfortable 
entering temperature. 

g. A damper should be placed in the supply duct. This can 
be used to reduce the flow of air during cold days and 
to start the ventilation flowing in the correct direc¬ 
tion. If the ventilation flows in the wrong direction 
in a system in which the supply duct is close to the 
floor and the exhaust is close to the ceiling, then 

the fresh air distribution will be quite poor. The 
ventilation can be started in the correct direction by 
closing the damper when the occupants enter, leaving it 
closed until their presence has warmed the shelter air 
sufficiently to insure that the shelter air is warmer 
than the outside temperature. Then the damper can be 
opened and the flow should go in the proper direction. 
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h. Tiie opening to the exhaust duct should be at the ceiling 

level. The warmest air in the shelter is close to the 

ceiling and b\ locating t.i .haust at the ceiling level 

the warm air can be drawn off. Further, since the 
ceiling is the warmest location in the shelter, exhausting 
the air from there will assist in providing a larger 
thermal differential between inside and outside, and 
thus improve the ventilation rate. 
i 0 The more ventilation supply ducts used in a shelter, 
the better the fresh air distribution. The supply 
ducts should be evenly spaced throughout the shelter. 
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